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Abstract. We apply the method of transitionless quantum driving for time- 
dependent quantum systems to spin systems. For a given Hamiltonian, the driving 
Hamiltonian is constructed so that the adiabatic states of the original system obey 
the Schrodinger equation. For several typical systems such as the XY spin chain, 
the driving Hamiltonian is constructed explicitly. We discuss possible interesting 
situations when the driving Hamiltonian becomes time independent and when the 
' driving Hamiltonian is equivalent to the original one. For many-body systems, a 

crucial problem occurs at the quantum phase transition point where the energy gap 
between the ground and first excited states becomes zero. We discuss how the defect 
can be circumvented in the present method. 
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1. Introduction 

Understanding tlie dynamics of a quantum system is a fundamental problem in 
quantum mechanics and is important for practical applications. With recent advances 
in experimental techniques, we need to control the dynamics in a high-precision 
way [U E]. Examples are seen in Bose-Einstein condensates in optical lattices [3], 
quantum computations |4j and so on. These advances motivate us to design the optimal 
Hamiltonian rather than solving the problem under a given Hamiltonian. 

Theoretically, several acceleration methods have been discussed recently: quantum 
brachistochrone [5l E], fast-forward method [TJ El |9], transitionless quantum driving 
(TQD) ^10] and Lewis- Riesenfeld invariant-based engineering pT|[T2]. In these methods, 
we discuss how fast one can reach a desired state starting from an initial one. We 
construct a time- dependent Hamiltonian so that the time evolution is achieved efficiently. 

These methods have been applied to systems with few degrees of freedom such as 
the Landau-Zener two-level system and a particle in a harmonic oscillator potential. In 
the present paper, we treat systems with many degrees of freedom by using the method 
of the TQD [10]. For a given Hamiltonian, we construct the driving Hamiltonian so 
that the adiabatic states of the original Hamiltonian obey the Schrodinger equation. 
The advantage of using this method is in its simplicity. The generalization to many- 
body systems is a straightforward task. The method has been tested experimentally in 
Bose-Einstein condensates in optical lattices and is shown to be robust against control 
parameter variations jl3] . 

Although the essential point of the method can be seen in small systems, we want to 
discuss the potential usefulness of the method when it is applied to many-body systems. 
We can find a lot of interesting phenomena in such systems. The effects of quantum 
fluctuations play important roles there. Since the method of the TQD manipulates 
the same effect, we expect that the present analysis reveals the nature of the quantum 
effects from a different aspect. As there are many interesting phenomena and a lot 
of techniques have been established, it is worth applying the present method to spin 
systems. 

Practically, the method will be most useful when the flnal state is a nontrivial 
unknown one. Starting from the initial trivial state, we evolve the system to the 
nontrivial state by using the quantum fluctuations. Such idea is known as the quantum 
adiabatic calculations or, more generally, the quantum annealing [HI [15], [16] . The major 
problem of this method is that inflnite times are required to reach the flnal state. By 
using the method of the TQD, we expect that such problem is improved. We can also 
flnd what kind of quantum fluctuations we should use for the time evolution. 

When we consider many-body systems, one of the most interesting phenomena 
that cannot be seen in few-body ones is the quantum phase transition. It is known 
that the adiabatic approximation fails at the phase transition point where the energy 
gap between the flrst and ground states becomes zero. This problem is not improved in 
the present method since the driving Hamiltonian is divergent at that point, which can 
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easily be understood from the general formula [TU]. However, we discuss that this is not 
a disaster and we can possibly circumvent the problem by using quantum fluctuations. 

The paper is organized as follows. In the next section, we review the method briefly 
and discuss possible extensions. Then, we discuss a two-level system in section [3l a two- 
spin system in section H] and the one-dimensional XY model in section [5l Section [6] is 
devoted to summary. 



2. Transitionless quantum driving 

We treat a time-dependent quantum Hamiltonian Ho{t). To obtain the corresponding 
state under some initial condition, we solve the Schrodinger equation, which is usually a 
formidable task especially for many-body systems. The analysis becomes considerably 
easier when we consider the slow evolution of the Hamiltonian. In that case, the 
instantaneous eigen-energies and -states defined as 

Ho{t)\n{t)) = E,,{t)\n{t)) (l) 

determine the state of the system. Here, n is the index denoting each state. Then, if we 
start the time evolution from one of the eigenstate \n{to)) at t = tQ, the adiabatic state 
is given by 



exp (^-zj'^dt'E^t') - J'^dt'{n{t')\^\n{t'))^ \n{t)). (2) 

The first term in the exponential represents the counterpart of the time-evolution factor 
in the stationary states and the second one is a phase factor that generates the geometric 
Berry phase [T7]. The adiabatic state is a good approximation as the solution of the 
Schrodinger equation when the following adiabatic condition is satisfied: 

|(m(t)|^|n(t))| 

(Emit) - E4t)y ^ ' ^ ^ 

where m and n denote different energy levels. 

Berry found a formula of the driving Hamiltonian such that transitions to other 
states do not occur [10] . This means that the adiabatic state becomes the exact solution 
of the Schodinger equation 

tj^\Mt)) = Hit)\Mt)). (4) 

The new Hamiltonian H{t) is different from the original one Hoit). Applying the time 
derivative operator to the adiabatic state ([2]), we can easily obtain H{t) = Ho(t) + Hi{t) 
with 

i^i(t)= E(l-Ht))(m(t)|)4|rn(t))(m(t)|. (5) 

We note that the time-derivative operator acts only on the state just behind the operator. 
If we consider the time-evolution under this new Hamiltonian H{t) starting from the 
initial state \n{tQ)), the state remains the same eigenstate of the original Hamiltonian 



Transitionless quantum driving for spin systems 4 

Ho{t). We can consider the time evolution with arbitrary speed, which accomphshes 
our purpose to go beyond the adiabatic evolution. 

The original formula ([5]) by Berry was derived under the condition that there are 
no degeneracies of the eigenstates. We used the orthonormal relation {m(t)\n(t)) = Smn 
to derive the formula. It is possible to extend the formula to systems with degenerate 
states. The adiabatic state for such systems was discussed in p^. We specify the 
instantaneous eigenstates as 

Ho{t)\n,fi,t) = E^{t)\n,fi,t). (6) 

Each eigenstate is labeled by two indices n and fi. It has an eigenenergy En and the 
number of values of /i represents the degeneracy. For the degenerate states we can 
choose the appropriate basis such that the orthonormal relation 

(n, /i, t\n', n', t) = 6nn'Sf,^' (7) 

is satisfied. Then, the adiabatic state is written as 

\Mt)) = Eexp (-^ fdt'Enit')) c^;:\t)\n, ^^,t) . (8) 

The coefficients c^"-* are determined below. Inserting this expression to the Schrodinger 
equation and using the adiabatic approximation, we obtain 

jA\t) = -T.{^,H,t\j^\n,f,'M;\t). (9) 
This differential equation is solved as 



c 



^"nt) = Ef^2w^?^(^o), (10) 



where uj^^, (t) is a matrix element of the matrix 



U^''\t) = Texp (^-z^*rft'A(")(t')) . (11) 

The symbol T represents the time-ordered product and the nonabelian gauge field matrix 
A^^^ has its element as 



^4"l(t) = (n,/i,t|^|n,/i',t). (12) 



d 



MM 



'dt 



In this adiabatic state, it is straightforward to apply the method of the TQD. After a 
similar calculation as the original case, we obtain 

Hi{t) = E -J2\^^'^^'t){n,u,t\^ i^\n,n,t){n,fi,t\. (13) 

We note that the formula (fT3l) is applied when the degeneracies are maintained 
throughout the time evolution. The method does not work well when the level crossing 
occurs. This can be understood from the alternative expression of the formula ([5]) as 

^.W- E l'W)^f|5f (»W|. (14) 
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We see that the driver Hamiltonian diverges at the level crossing point. This problem is 
crucial when we apply the method to the many-body systems where the level crossing 
occurs constantly. Even if we restrict ourselves to the time evolution of the ground 
state, the problem occurs at the quantum phase transition point where the energy gap 
between the ground and first excited states goes to zero. 

Before investigating this problem further, we discuss a possible extension of the 
original method. The Berry's formula is applied for any instantaneous energy levels. 
However, for many-body systems, we usually are interested in the ground state only. 
Then, if operators give nothing when they are applied to the ground state, they can be 
neglected from the beginning. If we only consider some specific energy level n, we can 
use the state-dependent driver Hamiltonian 

Ht\t) = (1 - |n(t))(n(t)|).||n(t))(n(t)| + (h.c), (15) 

where (h.c.) denotes the Hermite conjugate of the first term. Thus, we have an 
arbitrariness in choosing the driver Hamiltonian. This is a useful property for practical 
applications. 



3. Two-level system 



We can see how the method works well in a simple two level Hamiltonian 

Ho{t) = h('\t)-&, (16) 

where h^'^^t) is a magnetic field vector in three dimensional space and a are Pauli 
matrices. This case was considered in |T0]. Applying the general formula to this 
Hamiltonian, we obtain the driving Hamiltonian H{t) = Ho{t) + Hi{t) with 

1 



Hiit) 



/i(°)(t) X /i(°)(t) 



(17) 



2|/i(o)(t)|2 

where ho(t) represents the time derivative of h^^\t). Thus, the magnetic field to apply 
is given by h{t) = /iW(t) + (/iW(t) x /iW(t))/2|/iW(t)|2. 

This simple example tells us an important fact about the nature of quantum 
fluctuations. Even if we confine the magnetic field in a plane, e.g., h^'^^t) = 
{hi{t) , h2{t) , 0) , the driver Hamiltonian gives a field in the z direction. The quantum 
fluctuation effects unavoidably produce all kinds of operators which act on states in 
the Hilbert space. In other words, we need a complex matrix to represent the driving 
Hamiltonian. 

As a simple example, we consider an oscillating field 
/ ho cos ut \ 



h^'\t) 



\ 



ho sin 



I 



Then, the driving field has the same form as the original one as 

^ fiQ cos ut \ 
h{t) = Hq sinut 



\ 



(19) 
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where 

Thus, by considering the time evolution under the field ( IT9l) . we can obtain the adiabatic 
state of the original field ( !T8|) . 

It is interesting to consider some special cases. The first one is when the driving 
Hamiltonian becomes time independent. The condition h = gives 

hoit) + l^^Mt) X Mt) = 0. (22) 

We see that ^o(^) is perpendicular to /lo(^) aiid ho{t). In the case of the oscillating field 
(fT8|) . the condition = gives 

2{hl + hi) = ujh3. (23) 

Then, the driving Hamiltonian has the static field 

/0\ 

h{t) = . (24) 

U/ 

Driving the system by the static Hamiltonian H, we can find the adiabatic state of the 
original Hamiltonian. We note that we must choose one of the eigenstates of Hoito) as 
the initial condition which is not the eigenstate of H. 

Another interesting point to be used in the following is the existence of the 
fixed point where the driving Hamiltonian is equivalent to the original Hamiltonian: 
H{t) = Ho{t). This case is of particular interest since the adiabatic state becomes the 
exact solution of the Schrodinger equation. However, in the case of the two level system, 
the fixed-point condition does not give an interesting situation. By using the condition 
h{t) = ho{t), we obtain hoit) = /(t)/io(to) where f{t) is an arbitrary function of time 
with the initial condition /(to) = 1- This means that the magnetic field points to the 
same direction throughout the time evolution. In this case, the corresponding system is 
reduced to a static one by the time reparametrization dt' = f{t)dt. In the following, we 
find in some cases that the fixed-point condition gives a nontrivial result. 



4. Two-spin system 

As a preliminary calculation to many-body systems, we consider a two-spin system. 
This is the simplest many-body system and we can learn typical problems which occur 
in general many-body systems. Furthermore, the model is directly applied to two-qubit 
systems. 

The Hamiltonian is given by 

Ho{t) = Ut)alal + Jy{t)alal + h{t){al + a^), (25) 
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where cr^'^'^ and (72'^'^ represent Pauh matrices for two kinds of spins. We write this 
Hamiltonian by using the z-basis | mi, 7712) = 1^2)2 where (7f\m)i = m|m)i and 



cr||m) = m\m)2 with m = ±1. Arranging the basis states as | + 



+ 1, —1) and 



1, +1), we obtain a two block form 
( 2h J^-Jy 



\ 



Jx Jy 






-2h 












Jx ~l~ Jy 





-1), 



(26) 



The matrix in each block is understood as the Hamiltonian of a two-level system and we 
can apply the formula in equation ( |T71) . The second block represents a system with the 
magnetic field only in one direction and no further transformation is required. Applying 
the formula ( IT71) to the first block, we obtain the driver Hamiltonian in the space of the 
first block 

- jy) - Kh ^Jyl ( 







(27) 



This representation is transformed to that in the original full space as 

1 h(^Jx Jy) h(^Jx Jy 
2' 



Hi{t) 



Jy)' 



4/i2 + (4 

To obtain a more useful form, we consider the unitary rotation 
H{t) = U{t) (H,{t) + H,{t)) U\t), 

where 



U{t) = exp 



(28) 



(29) 



(30) 



The rotation is performed around the z-axis and the angle 6 is determined so that the 
transformed Hamiltonian has the same form as the original one. The condition is given 
by 



tan 261 



1 /i( Jx. - Jy) - h{X - jy 



Jx Jy 4/1^ -|- (t/j. Jy)"^ 

Since the unitary rotation has a time dependence, the transformed state 

\m) = uit)\m) 

obeys the Schrodinger equation with the Hamiltonian 

m = m + ^-^K + 02) 

= jx{t)ala^2 + jy{t)alal + h{t){al + a'^), 



where 



jx{t) + jy{t) = Jx{t) + Jy{t), 
(jx{t) - jy{t)) COS20(t) = J,(t) - Jy{t), 



h{t) = h{t) + 



(31) 
(32) 

(33) 
(34) 

(35) 
(36) 

(37) 
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Thus, the time-evolution of the adiabatic state of the original XY Hamiltonian f l25p is 
described by the same Hamiltonian with different coupling constants. If we consider the 
time evolution of a state for a given Hamiltonian H{t), the state must be the adiabatic 
state of the Hamiltonian Ho{t). 

We can also find in the present system the fixed point where Hi{t) becomes zero. 
This condition is achieved by the magnetic field 

Kt) = Ci-Ut) - Jy{t)), (38) 

where c is a constant. Although this condition is a trivial one in the first block of the 
matrix (|26|) . it is not in the whole space of the Hamiltonian. The system cannot be 
reduced to a static one by the time reparametrization. We expect that such a fixed 
point has a special meaning for the dynamical problem, which we have not understood 
yet. 



5. One-dimensional anisotropic XY model 

As the simplest many-body system showing nontrivial ground-state properties, we study 
the one- dimensional anisotropic XY model. We consider N kinds of spins and write the 
Hamiltonian as 

N 

Hoit) = - E {J^m^]+i + Jyit)^]^li - Ht)cr]) • (39) 
i=i 

Taking J^. = or Jy = 0, we can also study the transverse-field Ising model. In the 
static case, this model has a quantum phase transition at + Jy = h where the energy 
gap between the ground and first excited states goes to zero. 

This Hamiltonian is diagonalized by the method of fermionization. We define the 
Jordan- Wigner transformation fl9[ [20] 



^ (crj + = exp (^n a^a^ j a], 



(40) 



2a]aj - 1. (41) 
Here, the operators aj obey the fermionic anti-commutation relations 

Qia] + a^tti = 6ij , (42) 

ajOj + ajtti = 0. (43) 
We also introduce the Fourier transformation 

a,- = a„e*'?^'. (44) 

The momentum q takes a value between — vr and tt, and its discrete value depends on 
the parities of N and the number of fermions in each state. Since the value becomes 
continuous at the thermodynamic limit, we do not specify it explicitly here. Applying 
the transformation to the Hamiltonian, we obtain Ho{t) = J2q Ho{q,t) where 

Ho{q, t) = - 2 [( J^. + Jy) cos q- h] aid, 



q^-q 

i{Jx - Jy) sin q{d\a}-.q + Sgfl-g) - h. (45) 
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The Hamiltonian is block-diagonalized and each block is specified by the absolute value 
of q. In each block, the state is specified exphcitly by |0), ajalg|0), SglO) and a^q|0) 
where |0) is the vacuum of fermions defined by dq\0) = 0. The Hamiltonian is given by 
Hoiq,t) + Hoi-q,t) = H^^\q,t) © I2 + I2 © [-2(J^ + Jj,)cosg]l2 where 

^(1). ^^ ^ ( -2/^ - Jy) sing \ 

° ^ -2t{J, - Jy) sin q -4[( J, + Jy) cosq-h]-2h )' ^ ' 

We have a two-level system again and the driving Hamiltonian is constructed in the 
same way as the previous calculations. It is written as 

Ml). Mq,t) . / 1 \ 

Hi{q,t) = — - — smgl 1, (47) 



where 



T / .\ _ [(Jx + Jy) COSg - h]{X - jy) - {J^ - Jy)[{jx + jy) COS g - h] 



[{Jx + Jy) cosg - + ( Jj. - J^)2 sin^ q 

(48) 

In the full space of fermionic states, the Hamiltonian is written in terms of fermion 
operators as Hiit) = J2q Hi{q,t) with 

Hi (g, t) = :h^^tll sin g (alalg + a^qh^ . (49) 

The present purpose is accomplished if we represent this term by spin operators. 
However, the representation ( l49l) has a g-dependent coupling which gives nonlocal 
hopping of fermions. As a result, the corresponding spin representation has nonlocal 
and many-body interacting terms. Unfortunately, the method gives an unrealistic 
Hamiltonian. 

This problem does not arise if we restrict ourselves to a specific state. In most 
of applications, we are interested in the ground state of the Hamiltonian. The driving 
Hamiltonian is tuned so that the low-lying states are controlled properly. In the present 
case, the low-lying states are denoted by excitations of the modes at small g. The 
momentum g in the coupling Ji(g, t) is neglected and we have the driving Hamiltonian 
for the low-lying states 

^i(^) ~ -^H^i^^ (^^-'j + ^-9^9) ' (^0) 



where 



7 /,\ — ^"^^ ~^ '^y ^){Jx Jy) {Jx Jy){Jx + Jy h) 

^'^'^ - {Jx + Jy- h)- • ^''^ 

In this case, it is easy to go back to the spin representation and we obtain 

H,{t) = - Y: + (52) 

^ i=i 
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If we consider the unitary rotation as was done in the two-spin system, we can have as 
the total driving Hamiltonian 

TV _ 

m = - E (Jx.(t)^jaj^i + J.(t)^M+i - Ht)cr]) , (54) 
i=i 

where Jx(t), Jy{t) and h{t) are determined properly. 

We note that Ji{t) is divergent at the phase transition point + Jy = h as we 
can understand from the general formula (IHl). This problem can be circumvented if we 
choose the fixed point condition Ji{t) = 0. It gives the driving schedule 

h{t) = J,{t) + Jy{t) + C(j,.(t) - Jyit)), (55) 

where c is an arbitrary constant. In this schedule, the phase transition takes place at the 
isotropic point Jxif) = Jy(t). In the isotropic XY model, it is known that the singularity 
at the point = Jy = h is weak. It is interesting to see that the fixed-point condition 
chooses such a point properly to go beyond the phase boundary. 

6. Summary 

In summary, we have developed the method of the TQD and studied possible 
applications. We stress that there have not been so many applications of this method 
so far. In this work we enlarged the applicability of the method to various systems. To 
establish the usefulness of the present method it is necessary to consider more examples 
reflecting realistic experimental situations. 

When the method is applied to many-body systems, we have mainly three problems. 
First, we need to solve the eigenvalue problem for the original Hamiltonian which is 
generally a difficult task. Second, even if we can find the driving Hamiltonian, it has a 
very complicated form including nonlocal and many-body interaction terms. Third, the 
driving Hamiltonian diverges at the quantum phase transition point where the energy 
gap goes to zero. Concerning the first and second problems, we have shown that they 
are eased if we consider only a specific state such as the ground state. In that case, the 
driving Hamiltonian can be reduced to a simpler one. For the third problem, it is possible 
to avoid the phase transition by changing the schedule in a proper way. Although there 
is no general prescription, we have shown such an example in the one-dimensional XY 
model. 

The formula of the driving Hamiltonian indicates the existence of nontrivial 
conditions. Among them, the fixed-point condition is the most interesting one since the 
time dependence becomes essentially irrelevant. We expect that such special conditions 
become the key to understand the quantum dynamics. 

We stress that finding transitionless Hamiltonian is useful not only for direct 
applications but also for understanding the nature of quantum fiuctuations. In order 
to construct the driving Hamiltonian, we need to utilize the full operator space. For 
example, in a two-level system, we need a magnetic field not in a plane but in the 
full three dimensional space. In some applications such as the quantum annealing, we 
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usually control Ising spin systems by a transverse magnetic field in a single direction. 
The present method clearly indicates that this is not appropriate and we should use more 
different kinds of operators. We expect that the present method will be a guideline on 
how to construct the optimal Hamiltonian. 

Another interesting problem to be studied is the robustness of the transitionless 
Hamiltonian. It is important to know how much the time evolution is sensitive to 
the control parameters. This was discussed experimentally in |T3] and we need the 
theoretical ground to understand the result. This problem will be clarified in future 
studies. 
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